Instability of charge ordered states in doped antiferromagnets 
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We analyze the induced interactions between localized holes in weakly-doped Heisenberg antiferro- 
magnets due to the modification of the quantum zero point spin wave energy; i.e. the analogue of 
the Casimir effect. We show that this interaction is uniformly attractive and falls off as r~^'^'^^ in 
d dimensions. For "stripes", i.e parallel (d — f )-dimensional hypersurfaces of localized holes, the 
interaction energy per unit hyperarea is attractive and falls, generically, like r~'*. We argue that, in 
the absence of a long-range Coulomb repulsion between holes, this interaction leads to an instability 
of any charge-ordered state in the dilute doping limit. 



It is still not clear what happens when a dilute concen- 
tration of holes is introduced into a quantum Heisenberg 
antiferromagnet (AF). This is one of the central issues 
in the theory of correlated electronic systems, especially 
as it relates [0 to the high temperature cuprate super- 
conductors and related oxides. One class of proposals 
[Ij-Q holds that the result is a spatially inhomogeneous 
"charge ordered" ground state. Unfortunately, numerical 
analysis of the stability of such states is often inconclu- 
sive because the typical energy differences between states 
are small, and the Goldstone modes (spin-waves) produce 
finite size effects which decrease slowly with system size. 

The goal of this paper is to investigate the thermody- 
namic stability of charge-ordered states in short-range 
AF spin models in the dilute doping limit. We cal- 
culate the induced interaction between well-separated 
clusters of holes due to their modification of the spin 
wave spectrum, and find that it is uniformly attractive. 
Specifically, the asymptotic long distance (r ^ oo) in- 
teraction between two hole clusters (see Fig. |l|) is of 
the form £ ~ —JSr~'^'^'^^ for a d-dimensional, spin-S* 
Heisenberg quantum AF with exchange coupling J (see 
Eq. ^). For extended clusters of holes, under generic cir- 
cumstances, the dependence on separation (but not the 
absolute magnitude) of the interaction energy can be re- 
liably estimated by summing the pairwise hole-hole in- 
teraction over all pairs; for example, the interaction per 
unit hyperarea between parallel walls of localized holes 
(i.e. codimcnsion 1 hypersurfaces which, with the case 
of d = 2 in mind, we refer to as "stripes") falls with their 
separation as 8 ^—JS r~'^. 

As a consequence of this attraction, in the absence of a 
long-range Coulomb repulsion between holes, all charge- 
ordered states with sufficiently small hole concentration 
are unstable to phase separation, although it is possible 
that there exist non-vanishing hole densities for which 
charge ordered states are stable Remarkably, we 

find that the correct asymptotic form of the induced in- 
teractions cannot be obtained in any finite order of naive 
perturbation theory, because of the singular effect of a 
marginally bound (zero-energy) spin-wave state associ- 
ated with the fact that doped holes actually change the 
Hilbert space, by changing the number of spins. How- 



ever, a simple modified perturbation theory can be con- 
structed which qualitatively reproduces the exact spin- 
wave results in all cases we have tested. 
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FIG. 1. Clusters of localized holes are modeled by remov- 
ing the bonds connecting the spins to the rest of the system. 
The additional bonds inside the disconnected areas do not 
contribute to the interaction between hole clusters. 
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FIG. 2. A stripe of holes forming an anti-phase domain 
wall is modeled as a line of weak bonds in a perfect AF. 

The Model: The number of spin degrees of freedom 
changes with doping and, therefore, the Hilbert spaces 
appropriate to the doped and undoped AF are different. 
Therefore, for mathematical convenience it is preferable 
to use a model with a spin S operator on each site, and 
to treat any system with localized holes as a limiting 
case in which the coupling between a set of "impurity" 
sites and its neighbors goes to zero. For static holes, 
this is all there is to the model. However, so long as the 
holes are localized, either by an impurity potential or by 
a self-consistent field (as in Hartree-Fock solutions ||]) 
the effect of hole-hopping can be treated by including a 
larger set of modified exchange interactions in the neigh- 
borhood of each hole. Thus, the spin Hamiltonian of the 
doped system differs from that of the pure AF only in 
the strength of some exchange couplings: 



(1) 



where is the Hamiltonian for the perfect antiferro- 
magnet, which, for concreteness, we take to have only 
nearest-neighbor interactions on a hypercubic lattice. 
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and the perturbation Hamiltonian H' specifies a set of 
pairwise exchange interactions such that, in the limit 
A = 1, a spin near which the hole is localized is discon- 
nected from the rest of the system. Clearly, the interac- 
tion energy between hole clusters is obtained correctly in 
this limit, although the cluster self-energy could depend 
on the interactions between the fictitious, disconnected 
spins. We start with the simplest case of strongly local- 
ized holes, for which H' consists of the sum over the ex- 
change interactions connecting the designated hole sites 
in a cluster to the nearest-neighbor sites surrounding it. 
One exceptional geometry which we treat differently is 
a stripe which is simultaneously an anti-phase domain 
wall in the AF order. Such a stripe can be treated as 
shown in Fig. ^, effectively as a wall of bonds with altered 
exchange coupling, so that we work in the proper Hilbert 
space from the beginning; here naive perturbation theory 
in powers of A yields qualitatively correct generic results. 

Casimir Energy of Stripes: We begin by consider- 
ing the interaction energy between two stripes of static 
holes; this turns out to be the simplest problem because 
for A = 1 the region between the two stripes is cut off 
from its surroundings, and, according to a very general 
argument originally due to Casimir the interaction 
energy must fall off as £ ~ —r~'^. To be more explicit, 
if we take the stripes to be perfectly reflecting, the per- 
pendicular component of the spin-wave's momentum in 
the region of width r between the stripes is discretely 
quantized, which modifies its vacuum energy. Then the 
distance-dependent part of the energy per unit hyperarea 
can be expressed as the difference 



i2n) 



d-l 



rE^ 



where is the spin- wave energy, and Eoo is the vacuum 
energy density of the infinite lattice. The summation, 
performed with the Poisson formula, gives 
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which is proportional to the spin-wave velocity, c {fiiij^ ~ 
c\k\; in linear spin-wave approximation c = JS/Vd,) 
since only small momenta contribute to the large-dis- 
tance asymptotics of the interaction energy. This result 
is universal: it depends on the number of acoustic modes 
and their speed but not on the form of the spectrum at 
high energies or the specific boundary conditions. 

Linear Spin- Wave (LSW) Theory: In more gen- 
eral circumstances we need to compute the interactions 
from a more microscopic approach. For this purpose, we 



adopt LSW theory which is quantitatively accurate for 
large S, but which we expect to be a reliable method for 
extracting the long distance physics for d > 2 even for 
S = 1/2, since already in d = 2 AF order is very robust 
In order to make contact with the perturbative re- 
sults discussed below, we calculate the LSW correction 
to the ground-state energy in the presence of localized 
stripes or holes at arbitrary A. The exact ground-state 
energy can be calculated as a coupling constant integral 
of the expectation value of the perturbation Hamiltonian 
using the Feynman-Hellman formula, or, in the LSW ap- 
proximation, expressed as the functional determinant 
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The latter can be rewritten as a coupling-constant in- 
tegral of the diagonal part of the exact LSW Green's 
function (GF) as we shall see below. 

LSW theory for stripes: As a first application of LSW 
theory, we recompute the interaction energy between two 
stripes oriented along the y axis in d = 2. In an obvious 
mixed representation, labelled by the conserved wave vec- 
tor ky parallel to the stripes and the lattice position x 
in the direction perpendicular to the stripes, the Dyson 
equation for the GF reduces to a finite sum. 
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where i labels the vertical lines of sites connected by the 
weak bond representing the effect of the stripes, and the 
unperturbed GFs are given explicitly by 
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, (6) 



where the operators = ^r+^r+x, s,v = &J+^r+x are defined 
on the bonds (r, r-f-x) in terms of the Holstcin-Primakoff 
boson operators 5r, h\. As for the problem of a quantum 
particle in the presence of a finite number of point scat- 
terers, the solution of Eq. (|^) involves the inversion of 
only a finite matrix, and the functional determinant (^) 
can be rewritten in the convenient form. 
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with integrations over the coupling constant A, the zero- 
temperature Matsubara frequency lj, and the conserved 
momentum ky along the stripes, along with the finite 
summation over the stripe index i. 

The matrix elements of the pure crystal GF in the 
mixed representation are given explicitly by the integrals 
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where 7^ = [cos(fc^) +cos(fcj^)]/2 and = 1^7^- For the 
case of only two anti-phase stripes separated by the dis- 
tance r= 1x2 — xi I, the couphng-constant integration (|^) 
can be performed analyticaUy, with the result 
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The matrix elements of the bare spin-wave GF are 
small at r — > 00, so for all X<1, an asymptotic expression 
for the interaction energy can be obtained by expanding 
Eq. (p|) in powers of J>, D^', to leading order we obtain 
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The nai've perturbative results for this energy can be ob- 
tained by further expanding this expression to second 
order in A: 



£ 



(r)=- 



JSX^ 



153 



2048 7rV2. 



+ 0{r~^). (11) 



At A=l, however, the expression ( |lO| ) diverges, because 
the denominator of the argument of the logarithm in 
Eq. (^ has a zero at the point to — ky — O. It is this zero, 
identified as the zero-energy spin wave state bound to 
each stripe, which is ultimately responsible for the modi- 
fication of the asymptotic form of the interaction energy. 
As discussed later, the existence of such a state is re- 
lated to the change in the structure of the Hilbert space 
as the clusters become isolated at this value of A; such 
a state exists near a hole cluster of any geometry, and 
it cannot be eliminated by corrections due to spin-wave 
interactions or to the holes' mobility. In our calculation, 
we account for this zero-energy state by solving the scat- 
tering problem near each stripe or hole cluster exactly. 

The correct asymptotic behavior at A = 1 can be ob- 
tained by re-evaluating the interaction energy, starting 
with the complete expression (^). Using the long-distance 
asymptotics of the components of the GF, we obtain 
to leading order in 1/r 
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which leads to the universal Casimir result (||) evaluated 
at d = 2 and c = JS/ We have also investigated the 
crossover from the perturbative expressions (11), ( |lO| ) at 
small values of A to the universal form (H) at A = 1 by 
integrating the exact LSW energy (H) numerically. The 
results are shown in Fig. |3[ along with the corresponding 
asymptotic expressions. 

LSW Theory for Isolated Holes: Zero-energy spin- wave 
bound states and the associated divergence of the per- 
turbation series happen not only for extended objects 
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FIG. 3. The LSW Casimir energy of two stripes in units of 
JS'^{\^/S) from Eq. (9). The A = line indicates the pertur- 
bative result. The solid lines show the calculated asymptotic 
behavior at r — > 00 at intermediate values of A. 



like the stripes we just considered, but also for solitary 
holes or bigger hole clusters. As a second example, we 
consider explicitly the interactions between two isolated 
holes. The Hamiltonian of a single hole can be expressed 
as a sum of terms with different symmetries with respect 
to the point group and we find that in arbitrary di- 
mension, a zero-energy, totally symmetric bound-state 
emerges at the same critical value A = 1, so that the 
corresponding spin-wave scattering amplitude diverges at 
small frequencies. Thus, while for any A< 1, the interac- 
tion between two holes has the same asymptotic behavior 
as the leading order perturbative expression. 
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for A = 1 the asymptotic form of the interaction is 
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which is valid both for holes at odd separations (different 
sublattices, holes with opposite spins) and holes at even 
separations (same sublattice, holes with the same spin). 

Marginally -Bound States and the Proper Way to Do 
Perturbation Theory: We return now to the issue of the 
failure of naive perturbation expansion in powers of A. 
For a hole cluster of arbitrary geometry the Dyson equa- 
tion [Eq. (^ for the case of stripes] gives bound states 
when 
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where the two-component GF (g) is calculated for each 
pair of bonds present in the perturbation Hamiltonian 
H' as specified in Eq. (|^). At A=l, where the artificially 
introduced spins become disconnected, they acquire the 
freedom to rotate with respect to the rest of the system. 
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which reveals itself as a zero-energy spin-wave state lo- 
calized on these spins. Because of the continuity of the 
GF in Eq. (p^, this implies the existence of a soft spin- 
wave mode coupled to the holes, which is responsible for 
the divergence of the spin- wave scattering amplitudes off 
the hole cluster at small frequencies, and also provides 
the singularity required to modify the Casimir interac- 
tion energy obtained perturbatively. 

We checked this argument by computing the eigenval- 
ues of the matrix -y (w = 0) for several hole cluster 
geometries, including those shown in Fig. |1|. In agreement 
with the prediction, there are exactly two zero eigenval- 
ues for the clusters with only one disconnected group of 
spins, one for each spin wave branch, and there are four 
such eigenvalues for the geometry in the second example 
in Fig. 0, where the two spins of the cluster are discon- 
nected from each other. All other eigenvalues are pos- 
itive, implying that the perturbation expansion is free 
of singularities for A < 1. Within the same approach, 
we have also considered the effect of virtual hops of less 
strongly localized holes, which induce additional AF ex- 
change couplings between the spins in the neighborhood 
of the hole. Once again, the zero-energy spin-wave bound 
state persists regardless of the strengths of these addi- 
tional interactions. 

Because the existence of a zero-energy spin wave bound 
state within an isolated cluster representing localized 
holes follows from symmetry arguments, it is a very ro- 
bust feature of the spectrum. However, the degeneracy 
of the artificially inserted spins can be lifted by a local 
magnetic field term, added to the perturbation Hamilto- 
nian H' (obviously, this does not modify the physics of 
the rest of the system). It turns out that this term not 
only removes the divergence of the perturbation series, 
but also contributes an energy £ r^r^'^'^+^ to the interac- 
tion between the hole clusters already in second order of 
the perturbation expansion, which now gives the correct 
qualitative result. The improved perturbation expansion 
can be used to prove that the Casimir interaction be- 
tween hole clusters indeed falls off as £~r~^'^+^ even in 
the presence of spin-wave interactions. 

The only case where the non-physical spins need not be 
introduced is the stripe serving as an anti-phase domain 
wall, modeled as shown in Fig. ^ Although perturba- 
tion theory here also breaks down at A = 1 because of 
the soft spin-wave mode bound to the stripe, this mode 
is related to the freedom to rotate the AF magnetiza- 
tion of the two parts, which become disconnected at this 
value of A. Such symmetry is easily destroyed by the 
holes' mobility or by spin-wave interaction corrections, 
and, therefore, the robust asymptotic form of the inter- 
action energy between the antiphase stripes is the one 
given by the perturbation expansion (11). Note, how- 
ever, that the general asymptotic behavior £ ^~J' Sr^'^ 
is restored if next-nearest neighbor interactions are in- 
cluded in the Hamiltonian (0). 



Implications: We have considered Casimir interac- 
tions between well separated hole clusters in AFs. For 
hole clusters or stripes in a uniform AF, this energy is 
uniformly attractive and generally falls off with distance 
as £~r^^''+^ and £^r^'^ respectively. The interaction 
is quantitatively weak; for two holes in the S — 1/2 AF 
the interaction between next-nearest neighbor holes is 
less than 10^^ J. However, because the interaction falls 
slowly with distance, it is important for an analysis of 
the stability of static charge-ordered structures in sys- 
tems lacking long-range Coulomb repulsion. It has been 
conjectured [^jlo| that phase separation is a ubiquitous 
feature of lightly doped antiferromagnets, and that con- 
sequently there is always a first-order transition separat- 
ing the undoped and doped states. Evidence in support 
|ll| and in conflict |^,|l2| with this conjecture has been 
obtained from numerical studies of small-size systems. 
Phase separation has been shown to occur [^,0 in the 
large d limit of the Hubbard and t-J models, and in the 
mean-field spiral states of the large N t-J model 
The present results offer strong additional support for the 
validity of this conjecture. Specifically, we claim that be- 
cause of this Casimir-like interaction, any static ordered 
state of neutral holes will be thermodynamically unstable 
with respect to phase separation at small enough doping. 
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